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Abstract In order to find out whether empty singular boundaries can arise in 
higher dimensional Gravity, we study the solution of Einstein's equations con- 
sisting in a (N + 2) -dimensional static and hyperplane symmetric perfect fluid 
satisfying the equation of state p = r/p, being p an arbitrary constant and N > 2. 

We show that this spacetime has some weird properties. In particular, in the 
case rj > —1, it has an empty (without matter) repulsive singular boundary. 

We also study the behavior of geodesies and the Cauchy problem for the prop- 
agation of massless scalar field in this spacetime. For r) > 1, we find that only 
vertical null geodesies touch the boundary and bounce, and all of them start and 
finish at z = oo; whereas non- vertical null as well as all time- like ones are bounded 
between two planes determined by initial conditions. We obtain that the Cauchy 
problem for the propagation of a massless scalar field is well-posed and waves are 
completely reflected at the singularity, if we only demand the waves to have finite 
energy, although no boundary condition is required. 

Keywords Gravitation • singularities 
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1 Introduction 

We have recently pointed out that solutions of the Einstein field equations with 
empty singular boundaries can exist in 4 dimensions 1 , 2, 3 , 4J . These singularities 
are not the sources of the fields, but rather, they arise due to the attraction of 
distant matter. 

In such spacetimes, at the boundary, some components of the Weyl tensor di- 
verge while the energy-momentum tensor vanishes. It turns out that only orthog- 
onal null geodesies just touch the singularity and bounce, whereas non-orthogonal 
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null ones as well as all time- like ones bounce before getting to it. Furthermore, 
we have shown that the Cauchy problem for wave propagation is well-posed, if we 
only demand the waves to have finite energy, although no boundary condition is 
required [5]. And waves are completely reflected at the singularity. Due to these 
properties we also call this kind of singularities white wall^\. 

Needless to say, the possibility that space-time may have more than four dimen- 
sions is now a standard assumption in High Energy Physics. Indeed, most theories 
attempting to unify all the fundamental interactions require a higher dimensional 
spacetime. 

However, it is well known that Einstein's Gravity gets weaker as the number 
of spacetime dimensions increases. For instance, stable bound Keplerian orbits 
do not exist in spacetimes of more than 4 dimensions [6|, or also, the degree of 
compactification of spherical stars diminishes as the dimensionality of spacetime 
raises [7]. 

Thus, it would be worthwhile finding out whether empty singular boundaries 
can arise in higher dimensional Gravity, and the aim of this paper is to explore 
this issue. 

In this work, we show that such singularities may occur in spacetimes of TV + 2 
dimensions (TV > 2) by finding an explicit solution. 

The solution to be described is a (TV + 2)-dimensional static and hyperplane 
symmetric spacetime filled with a perfect fluid satisfying the equation of state 
p = rjp, being r\ an arbitrary constant. 

In section [2] we find the solution and describe its main features depending on 
the value of r\ and TV. In particular we show that, for any TV > 2, this spacetime 
has an empty singular boundary if r) > — 1. We also discuss the vacuum and dust 
limits. 

In section [3] we study the properties of null and time-like geodesies for the 
most physically interesting case r\ > 1. 

In section [H we discuss well-poseness of Cauchy problem for the propagation 
of a massless field in this case. 

Throughout this paper, we adopt the convention in which the spacetime metric 
has signature (— + ■ ■ ■ + +), the system of units in which the speed of light c = 1, 
and G is the (TV + 2)-dimensional Newton's gravitational constant. 



2 The solution 

Let us consider a solution of Einstein's field equations corresponding to a static 
and hyperplane symmetric distribution of matter. That is, they must be invariant 
under n translations and under rotations in TV(TV — l)/2 planes. The matter we 
will consider is a perfect fluid satisfying the equation of state 



where rj is an arbitrary constant. Then the components of the stress-energy tensor 
read 



P = 



VP, 



(1) 



Tab =P(V+ 1 ) u aUb +Pg a b, 



(2) 



A reviewer suggested for them the name singular mirrors instead white walls. 
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where u a is the velocity of fluid elements. 

Due to the required symmetry, by solving Killing's equations it can be readily 
seen that coordinates can chosen such that the line element can be written as 

ds 2 = -G{zf dt 2 + e 2V ^ ((dx 1 ) 2 + ■■■ + (dx N ) 2 ) + dz 2 . (3) 
The non identically vanishing components of the Einstein tensor are 

G tt = -NG 2 (v"+ ^±1 



G NN = e 2V (9l + [N- 1) J V' + (N- 1)V" + 



V' 2 




V 2 ) 


,(5) 


M 


,(6) 



G„ = NV ^ + — 

where a prime (') denotes differentiation with respect to z. 

On the other hand, since the fluid must be static, u a = {—0, 0, . . . , 0), so 

-.2 2V 2V 



Tab = diag [pG ,pe ,...,pe ,p\ , (7) 

where p and p depend only on z. Note that due to the imposed symmetry, rotation, 
expansion and shear of the fluid velocity vanish, while the acceleration vector field 

is u a = (p,...,Q,g'/g). 

Thus, Einstein's equations, i.e., G a b — 4^G JV W 1 T^bEl become 

V» + "±±V» = -£^ P , (8) 

% + tN-l)*V + (N- D V + = ^P, 0) 

*V' + ^V» = glL P . (10) 

Moreover, V a T ab = yields 

p' = -(p + p )<L. (11) 

Of course, due to Bianchi's identities equations (f8)l- (|lip are not independent. 

Some comments are in order. First, it is interesting to see that the (N + 2)- 
dimensional Raychaudhuri equation 

+ 6 2 - V a u a + a ab a ab - co ab co ab + AtvG (p + jj±jp\ = , (12) 

becomes in this case 

V a u a = G ^+N G -V' = A, G L+?L±± P \ , (13) 



2 See appendix 1X1 

3 We have set the normalization constant such that, in the Newtonian limit, they lead to 



;V 2 (l + 3tt) = 4nGp. 
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which is equivalent to the linear combination — (N— l')(f8)l + (f9))+A r (|10|) of equations 

©-CG3I. 

Second, it readily follows from (fg))- (|lip that, for arbitrary constants A, B, X ^ 
and zo, field equations are invariant under the transformation z — > Xz + zq (i.e., 
^-translations, z-dilatations and mirror reflections across any "plane"), Q — > AQ, 
V^-V + B, p -> X 2 p and p -> X 2 p. Thus, if {G{z),V(z),p(z),p{z)} is a solution of 
©-(HO]) {AG(Xz + z ),V(Xz + zo) + B, X 2 p(Xz + z ), X 2 p(Xz + 2 )} is another one. 

Third, regarding the mirror symmetry, independently of the equation of state 
and dimension of spacetime, a solution cannot have a "plane" of symmetry in a 
region where p(z) ^ and p{z) > 0. In order to see this, let us assume that z = z s 
is that "plane" , then it must hold that Q' = V' = p' = p' = at z s , and so from 
(jTDJ) we get that also p(z s ) = 0. Now, by differentiating (fTTj) and using ([Sj) and ([9]), 
we obtain p"(z s ) = —AnGp 2 < 0. By using the equation of state we can go further. 
Indeed, it follows from (fTj).(f8))- (|10|) that the four functions Q, V, p and p, and all 
their derivatives of any order must vanish at z s . 

We now proceed to find the solution. By using the equation of state ([1]), we 
readily find from that 

p = po g-(v+D j (14) 

where po is an arbitrary constant. Now, this last expression and the change of 
variables 



f g{z')-^ -dz 1 ', (15) 

J Zn 



brings © and (10]) to 



rj+1 G$ N + l 2 4ttG 

--2--g v e + v ee + -2- v e=-N=iiP° ( 16 ) 

and 

Gj N-l 2 AtvG , . 

-gV 6 + ^-Vs =J^JP0- (17) 

Hence, we find that V satisfies 

r,(N - 1) + 3iV + 1 a 2tvG 
% + 4 ^ + ]^3T W " 1) Po = °. ( 18 ) 

and, by setting 

(q-l)(^(JV-l)+3JV + l) 

a= 2(ATl) " GP0 ' (19) 

we can write its general solution as 

V(£) = In (Ci sin(^7i£ + C* 2 )) " (N - 1 » +3JV+1 , (20) 

where C\ and C2 are arbitrary constants. 
Now, we can write (|17p as 

di ln{G) -N^TV^^- V t' (21) 

which, by using (|20[) . can be readily integrated, giving the general solution for (?(£) 

_ 2(AT-1) _ 2_ 

5(0= (C3sin(V3C + Ca)) " ( ™- 1,+3N+1 (cos(VTi£ + C 2 )) , (22) 
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where C3 is a new arbitrary constant. 

Now, we define a new positive coordinate u such that 

sin 2 (^c7C + C , 2 ) 
u = o 1 ( 23 ) 

and, without losing generality, we set C\ = C3 = l/-\ffi, being j3 is an arbitrary 
constant. And, taking into account that 

du du d£ 1 4a 1 . ,1 __2±I 

i = iT u 9 2 ' (24) 



we can write the the metric ([3]) in terms of - 



2(JV-1) 2 „ 4 / , 9 0\ 

V (n-i)+3n+i (1 _ p u ) v -i dt + u^ N ~ 1 )+ 3N + 1 Udx ) H h (da; ) J 



ds 

4- 

4a 



ft 2^(W-l) + 4JV 2,7 , 

+ -T- u '(»-i)+3»ti (l-/3ti) du 



/ 1^ ^ «^ J<1//? if/J>0 

—00 < t < 00, —00 < x < 00, . . . , —00 < x < 00, < it < . (25) 

|<oo if/3<0 V ^ 

For the case N = 2, this solution was found in [8], see also [9l ll0lfTT| : and a more 
detailed study of its properties was done in [2]. Note that, for a Lorentz manifold, 
a and /3 must have the same sign and none of them can be zero or infinity (see 
also ([23])), therefore < a/3 < 00. On the other hand, it follows from (|23[) that, if 
a > the coordinate u is bounded (0 < u < 1 //3) , whereas if a < it is unbounded 
(0 < u < 00). 

Note that, by rescaling the coordinates, the parameter a might be eliminated. 
However, we prefer to keep it and think of it as a scale with dimensions of 
(length)^ 2 . In this way, the parameter /3 and the coordinate u turn out to be 
dimensionless. In this sense, for each N > 1, the solution essentially depends on 
only two dimensionless parameters /3 and 77. 

Even though in this work we are mainly interested in the case 77 > 1, it is worth 
mentioning that t)25|) has a much wider range of validity. Indeed, for each N > 1, 
from the mathematical point of view the solution given in (|25p makes sense for 



arbitrary finite values of the parameters, provided a/3 > 0, r) 1 and r\ ^ — 3 ff^ ■ 
Even more, for these excluded values of 77 the solution can be found as limit cases 
of (|25[) . Nevertheless, for the sake of completeness, we shall give here a briefly 
description of the main features of the solution for arbitrary values of 77. 

In what follows, let us assume that < a/3 < 00, |?7| < 00, 77 ^ 1 and 77 ^= — i^rr ■ 
In view of (114[) , (|19[) and (|25|) , the fluid density may be written as 

2n(N — 1) a 17+1 
^ )= ^ )= (,-l)(/l) + Ul) tf «^-^ (1- (26) 

Hence, the sign a must to be chosen so that 077(77 — l)(ri(N — 1) + 37V + l) has 
the same sign as the fluid density, and this choice sets the sign of /3 as well. 

Therefore, for p > 0, the parameter a is positive for the cases 77 > 1 or > 
77 > — 3 $^i , but negative for 1 > 77 > or — 3 ff-i > V- Thus the range of the 
coordinate u is (0 < u < 1//3) for the former cases and (0 < u < 00) for the latter 
ones. On the other hand, it follows from (|26|) that p(u) and p(u) are C°° functions 
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of u in the interior of the range. The pressure is positive for 77 > and negative 
for > 77, and the dominant energy condition \p\/p < 1 holds only for 77 > 1 or 
-1 > 77. 

For 77 = —1, which is a case with a > 0, it follows from (|26[) that p = —p = 
2(n+i)ttG • Hence, in this case, (|25[) turns out to be a vacuum solution with cos- 
mological constant A = (jy+i) Q > * ms ^ s * ne higher-dimensional generalization of 
that found in [12] , 

For 77 ^= — 1, it follows from (|26[) that at the endpoints of the range of u, 
depending on the sign of the exponents of u and (1 — /3u), p(u) either vanishes 
or diverges. At u = the exponent is ^(jv-V/^jv+i ' so ^ vanishes for the cases 
77 > —1 or — 3 x^i > 77, but diverges for — 1 > 77 > — 3 ^ ±^ ■ Thus, this property 
together with the sign of a define five regions where the solutions with positive p 
have very different features. Let us enumerate these regions in the following way, 
(I): 77 > 1, (II): 1 > 77 > 0, (III): > 77 > -1, (IV): -1 > 77 > and (V): 

— > 77. In Table[T]we list the main properties of the solution for the different 

regions. 

Similarly it follows from (|26[) that, at the other endpoint (u = 1//3 or u = 00, 
depending on the sign of a), p(u) vanishes for the cases (I), (II) and (IV), but 
diverges for the (III) and (V) ones (see Table [1}. 

Note that the proper distance coordinate z used in ([3]), is implicitly given in 
terms of u by 



±/ \/guu(u')du =±J f- / (m'^-d+sjv+i- 1 (i-p u 'y~~ 



du 



C + C' M *7(N-i)+3w+i (i + o(u)) asit^O, 

C + C'{l~l3uy^ (1 + 0(1 - Pu)) for 13 > as 1/(3, (27) 

C + C u i7(«-i)+3N+i v -i (l + 0(u -1 )) for f3 < as u 00, 



where C and C' are given constants for each of the cases. 

When u — y it follows from (|2T)l that, z -¥ z in cases (I), (II), (III) and (IV), 
and z — > ±00 in case (V). At the other endpoint, z — ► ±cxj in cases (I) and (II), 
and z — > z\ in cases (III), (IV) and (V). Therefore, spacetime is semi-infinite in the 
z-direction in cases (I), (II) and (V), while it is bounded between two boundaries 
in cases (III) and (IV). 

The solution may have a spacetime curvature singularity only at the endpoints 
of the range of u, since straightforward computation of the Kretschmann curvature 
invariant yields 



K = R 



abed 



r: 



abed 



4(W+i) . 

256 a 1 u .,<«-i)+3jv+i (1 _ p u yr,-i 



N 2 (N 2 -i)( v -iy 



/9 2 (V ~ I) 2 (v(N - 1) + 3iV + 1) 
+2N 2 (N- 1)(t? 2 - 1) ((N- 1)77 + N + 1) (3u + P(N, V ) p 2 u 2 ^ , (28) 
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Region 


I 


II 


III 


IV 


V 


V 


77 > 1 


1 > rj > 


> f? > -1 


1 > n > 3iv+1 


3JV+1 
JV-1 


a,$ 


+ 




+ 


+ 




range of u 


(0, 1//J) 


(0,oo) 


(0.1/0) 


(0, 1/j8) 


(0,co) 


P 


+ 


+ 


- 


- 


- 


\P\/P 


< 1 


> 1 


> 1 


< 1 


< 1 


u ->■ 0+ 


p -> 


p -> 


p -> 


p — > —00 


p -> 


u (1/j8)- 

u — ¥ CO 


p -> 


— 
p -»■ 


p — > —00 


p -> 


— 
p — > —00 


u -> 


-> zo 


-> zo 


-> zo 


-+ ^0 


-> ±00 


« -> 1/0 

M — > OO 


->■ ±00 


-> ±00 


-> Zl 


-> Zl 




w = 


singular 


singular 


singular 


singular 


regular 


u = l//3 

M = OO 


regular 


regular 


singular 


regular 


singular 



Table 1 Some properties of solutions with positive p according to the value of r\. 



where P(N, rj) is the polynomial 

P(JV, 77) = A 4 + 9A 3 + UN 2 + 7 A + 1 + 4(A - 1) (A 3 + 6A 2 + 5A + l) 77 

+2 (a 2 (A - 1) (37V + 8) + A 2 + AN + 3) r? 2 

+(A - 1) (A(A - 1) + 1) ((A - 1) 77 + 4(A + 1)) 77 s . (29) 

When u ->• it then follows from ([28| that, A' — > 00 in cases (I), (II), (III) and (IV), 
and K — > in case (V). At the other endpoint, K -> in cases (I), (II) and (IV), 
and z — > 00 in cases (III) and (V). Therefore, cases (I), (II), (IV) and (V) have 
one singularity, and case (III) is a spacetime trapped between two singularities. 

Note that at the singularity at u = 0, p,p — > in cases (I), (II) and (III), and 
these singularities are of the kind that, because of matter vanishes at them, we 
call empty singularities. 

In contrast, at the other three singularities, i.e., u = 00 in cases (II) and (V) 
and ii — O in case (IV), p and p diverge. 

It is easy to see that the singularities at u = are repulsive. Indeed, if we 
consider a test particle initially at rest near the singularity, we get from (|25[) that 

= -\g uu a»in|$tt|>o (30) 



d 2 u 

dT 2 
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if „ > -3JV+1 
11 V -* 77-1 • 

In our discussion, we have left aside the four values of 77 that are at the end- 
points of the regions, for the sake of readability and completeness, the solutions 
for these cases are postponed to appendix [Bl 

It follows from ([TJ, (]14[) and (|19[) that p and p are proportional to a, then 
we can get a vacuum solution from (|25[) if it limit a — >■ makes sense. But, we 
immediately see that this only occurs if also /3 — > keeping their ratio fixed. Thus, 
proceeding in this way and after the change of coordinates 

1 + HZ = U 1(K~ N 1) + 3N + 1 t (31) 

we get from (|25[) the vacuum solution 

da 2 = - (1 + k2)~ 2 (^Tt) dt 2 + (1 + Kz) 1 ^(^dx 1 ) 2 H h (dx™) 2 ) + dz 2 , 

— OO < t < OO, — 00 < x 1 < OO, . . . , — OO < X™ < OO, < 1 + KZ < oo , (32) 

where k is an arbitrary constant. This is the generalization to higher dimensions 
|13j of Taub's plane vacuum solution [13]. The new limit n — > 0, brings (|32p to 
Minkowski spacetime. 

Regarding the dust limit, it follows from (UJ that it corresponds to 77 — > 00. In 
this limit, (1251) becomes 



ds 2 = -dt 2 + (da; 1 ) 2 + • • ■ + (dx n ) 2 + vT 2 (1 - /3tj) -2 du 2 , (33) 



and the change of coordinate z = \J ^ In f 1 _^ M J shows that this is Minkowski 

spacetime. On the other hand, it follows from (|26[) that p vanishes as ?i -2 and p 
as t; -1 in this limit. Therefore, as expected, of course, there are not static dust 
solutions. 

Summarizing, we have studied the main properties of the solutions with posi- 
tive p as a function of the parameter 77, these properties are listed in Table [1] In 
particular, we have found that empty repulsive singularities appear in the cases 
with 77 > — 1. In a similar way, of course, we may analyze the less interesting cases 
with negative p, but we are not going to do it here. 

From now on, we restrict ourselves to the most physically interesting case I 
(77 > 1). We recall that, in this case, spacetime is semi-infinite and the range of the 
coordinate u is < u < 1//3. On the other hand, (|28|) shows that there a singularity 
at u = and K — > as u — > 1//3 (z — > 00) suggesting it is asymptotically flat at 
spatial infinite in the z direction. It follows from (|2(jp that pressure vanishes at 
u = and u = 1//3, it is positive elsewhere and it has a maximum at 

— ^-l){N-l) (34) 



2/3(t7(JV-1)+JV+1) ' 

Since p and p vanish at the singularity, we see that it is an empty singularity. 
Thus, the attraction of matter curves spacetime in such a way that a repulsive 
singularity arises in a place free of matter. 

In the following sections we shall study movement of particles and propagation 
of waves in this spacetime. 
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3 The geodesies 

In this section we study the geodesies in this spacetime, in the case r\ > 1. Since the 
components of the metric are independent of t, x 1 , . . . , and x N , the momentum 
covector components pt, pi, and p N are constant along the geodesies. For 
timelike geodesies, we set fi 2 = 1, E = —pt/m and pi = Pi/m, and we choose r to 
be the proper time; and for null ones, we set fj, 2 = 0, E = —pt and pi = Pi, and we 
choose r to be an affine parameter. So, we can write 

rjf 2(W-1) 2 _ 

= u ,(jv-i)+3jv+i (i _ j3 u )— e, (35) 

dr 

dc£^ 4 

— = u ^(iv-i)+3N+i g l<i<N, (36) 

dr 

du\ 2 4 a 2^(w-i)+4w 



— = — ui(«- 1 )+ 3N + 1 (l-Bu)"- 1 
dr J p 

x ^(n~i) 2 +3«+i (l-pu)^ 1 E 2 - p 2 - u ~i(«-i)+3»+i p 2 j j (37) 

where p 2 = p 2 H + p N 2 . 

The right hand side of (|37[) cannot be negative, then it must hold that 



-o II T,(JV-1) + 3JV4-1 ,.f _|_ ^ 

E 2 > V(u) := — M +P 2 . (38) 

u „(iV-l) + 3N + l (1 _ ^ 

Therefore, only vertical null geodesies touch the singularity at u = (z = zo) and 
reach u = 1//3 [z = oo). Whereas, as it is shown below, non- vertical null ones 
(p 2 > 0) as well as massive particles bounce before getting to the singularity and 
bounce again before reaching z = oo. 

For the former case, the geodesic equation can be integrated in closed form. 
Indeed, when ^i 2 = and p 2 = 0, we get from (f3"5]l and (f3"7|) 

— =±JJ±- ^(iv-D+aiv+i- 1 (i-^)- 1 (39) 

thus (see, for example, [15 ) we have 

/IT n(N-l) + SN + l 3W-2 

t-tn = A/— — — u >7(JV-l)+3iV + l 

r,/ 37V -2 . , , 3N-2 a \ , An . 

x 2Fl ( ,(/v-i) + 3iv + r 1; 1 + , ( iv-i) + 3iv + i ; Pu ) ■ {A V 

Note that 



|* — *o| oc < (41) 
I — ln(l — j3u) — > oo as u — > 1/(3 ■ 

Therefore, a freely falling photon bounces off the singularity at time to, and then 
continues its upward travel to u = 1//3 (z = oo). 

It follows from (|38|) that the movement of non-vertical [p 2 > 0) photons or 
massive particles is constrained to the region where E 2 > V(u). Since, p 2 and fj, 2 
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are not both vanishing in these cases, the function V(w) is a positive continuous 
function of u for < u < 1//3 and, since r/ > 1 and N > 2, V(u) — > +oo when u — > 
oru-> 1//3. 

Moreover, V(u) turns out to be a convex function in this interval. It can be 
seen as follows. We can write V(u) as linear combination, with non negative coef- 
ficients, of two functions of the sort f(u) = u~ a (l — /3u)~ b with a,b > 0. Now, a 
straightforward computation yields 

f"(u) = U ~ a ~ 2{1 a ~ P b U) ~ " 2 (ab+(l + a + b)(a-(a + b)p u y^ , (42) 

which is clearly positive in < u < 1//3. Then, V"(w) > in this interval, and so 
there is one and only one point uq where V(u) attains a positive minimum. 

Note that, for massive particles, the value of mo depends on p 2 besides the 
parameters of the solution (N, rj and /?); whereas for photons with p 2 > 0, mo is 
independent of p 2 . 

However, for any case and any value of E 2 > V(tto), we find two (and only two) 
turning points u min and u m ax where E 2 = V(u). 

In the case of massive particles which are not moving horizontally (/i 2 = 1 and 
p 2 = 0) , we immediately see that u = uq is a stable equilibrium "plane" . So a 
particle at rest in this place will remain at rest, and E 2 = V(uo)- On the other 
hand, if E 2 > V(«o) the particle bobs up and down between u min and u ma x- 

In the case of non- vertical geodesies (p 2 > 0), if E 2 > V(uq) timelike geodesies 
as well as null ones are wiggly worldlines where the coordinate u oscillates be- 
tween u m i n an d u m ax, while the horizontal movement is unbounded. Whereas, if 
E 2 = V(«o)j the movement is horizontal and the geodesies are straight lines in 
the submanifold u = uq. However, as mentioned above, the position of the equilib- 
rium "plane" uo where horizontal geodesies occur, is independent of p for massless 
particles, but depends on it for massive ones. Nevertheless, it is not difficult to 
see, that the worldlines of horizontally-moving massive particles approach to the 
massless ones in the hight energy limit (p 2 3> 1). 



4 Reflection of waves at the singularity 

After finding that all geodesies bounce off the singular boundary, let us turn our 
attention to propagation of waves. In particular, we shall explore the behavior of 
an incident wave pulse when it reaches the singularity. 

To this end, in this spacetime, we consider the propagation of a massless scalar 
field with Lagrangian density 

^ = -\v\V aV = -\g ah da^d bV . (43) 

As usual, we obtain the field equations by requiring that the action 

S = J ^{Va V ,f,gab)V\9\dtdzd N x (44) 

be stationary under arbitrary variations of the fields Sip in the interior of any 
compact region, but vanishing at its boundary. Thus, we have 



Higher-dimensional perfect fluids and empty singular boundaries 



11 



which, in our case, reads 



da (V\g\g ab d b <p) 

V a VV = — i = J -=Q. (45) 



By using (I25|) we can explicitly write down the wave equation. However, it suffices 
for our purposes to consider the behavior of waves propagating in a neighborhood 
of the singularity. Thus, taking into account that near the singularity /3u <C 1, we 
can therefore neglect the terms j3u in (]25f) . If, in addition, we make the change of 

JV-l , _ N 2JV 

coordinates t — > C 2 t, x % — > C x 1 and u — s> z = C "«i(»-w»ti ; where the 
constant C is given by 

C=( ^ d , (46) 

\I3( V (N -1) + 3N + 1) 



we get that near the singularity the metric (|25|) behaves as 

ds 2 w -(v^z)" 1 "^ (di 2 - dz 2 ) + (v^z)^ f (da; 1 ) 2 + ■ • ■ + (da;™) 2 ) . (47) 



The line element in the right hand side corresponds to the vacuum solution (|32[) . 

1 2 

this can be readily seen through the change of coordinate (-y/az)" — > (1 + kz) n + 1 . 

By using (|45|l and (147[) we get the wave equation satisfied by ip near the sin- 
gularity 

dttf = -d z (zd z <p) + — — — — —A N ip , (48) 

2 [y/az)^ N 

where is the usual Laplacian of the iV-dimensional "horizontal" space. 

The solutions of the wave equation (|48p have been studied in [5 . We briefly 
describe here the main results. We start with defining the underlying elliptic dif- 
ferential operator A and the Hilbert space H on which A is symmetric. 

Let Q = x (0, oo) and dfj, the Lebesgue measure on Q. Define, when <p £ 
Co°(n), the operator A by 

Ay = —8, (z d zlfi ) - 1 A<p. (49) 

Then consider the Hilbert space 

H := L 2 {f2,zdn) 

= {tp(x,z) : I \p(x, z)\ 2 z d/i < oo}. 
■In 

By construction the operator A is symmetric on H, with 

< A<p,r) > H = / f d z ipd z f] + — — — — —Vip ■ Vfj ) zd/j 
Jn \ {^z) L+ n 

for vp, 7? € Cg°(ft). 
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This leads to introducing the "energy space" 

£ = & 6 Hl oc {n) n H : b(ip, ip) < oo}, (50) 

where Hl oc (Q) is the usual local Sobolev space. It is straightforward to check that 
£, equipped with its natural norm 

llvlll : = 6 (v>> f) + IMIh, 

is a Hilbert space. This is the largest subspace of H on which the form b is finite 
everywhere. 

At the impossibility of provide a boundary condition at the singularity, our 
first question is whether A is essentially self-adjoint or not: as a result, it is not. 
Indeed, we probed that it holds [5] 

Theorem The operator A is not essentially self-adjoint. However, there exists only 
one self-adjoint extension of A whose domain D is included in the energy space £. 

However, we are only looking for those extensions with domain included in 
the energy space, because we are interested in waves having finite energy. When 
taking into account this restriction, we recover the uniqueness of the self-adjoint 
extension of A. The domain of this particular extension is 

D:={tpe£ :3C>0Vt)££, \b(<p,v)\ < C\\r)\\ H }. (51) 
Note that, there is no boundary condition attached to the definition of A. 

Now, coming back to the wave equation, we take suitable functions / and g on 
Q and consider the Cauchy problem 



dmp{t, x 1 , . . . , x , z) = —A<p(t, x 1 , . . . , x N , z), 
p{0,x\...,x N ,z) = f{x\...,x N ,z), (52) 
^(O,^ 1 , ...,x N ,z) = g(x\.. .,x N ,z). 

We found in [5] that, if / 6 £ and g S H, the problem (|52|) has a unique 
solution, and there exists a constant C > such that 

Vt>0 \Mt,-)\\ £ + \\dMt,-)\\H<C(\\f\\ £ + \\g\\ H ). (53) 

Furthermore, in this case, the energy 

Efat) := i / ( (cM 2 + (cM 2 + * |V^| 2 ) zd^ (54) 

is well-defined and conserved: 

Vt>0 E&,t) = ~(\\g\\ 2 H + b(f,f)y (55) 

This result shows that the Cauchy problem (|52|) is well posed without any 
boundary condition on tp. This does not necessarily mean, however, that if vanishes, 
or has no limit at all, at the boundary. Indeed, provided / and g are regular enough, 
ip does have a trace on dfl at each time t > 0, which is entirely determined by the 
Cauchy data (see [5], for details). 



Higher-dimensional perfect fluids and empty singular boundaries 



13 



In other words, if the energy of the initial field configuration is finite, it is 
conserved and, at any later time, the field is completely determined by the Cauchy 
data anywhere — even at the boundary. Therefore, no boundary condition should 
and can be provided. 

It is important to remark, that this result may be rigorously obtained by using 
Theorem 7.1 of [5], without making any approximation for ftu <g 1. If fact, by 
using (|25[) and (|45[) we can explicitly write down the exact wave equation. The 

relevant coefficient \f\g~\g uu , turns out to be J~^f- u(l — /3u), and its inverse is not 

integrable at the origen. 

Since the energy is conserved we see that waves are completely reflected at the 
empty singular boundary. 



5 Concluding remarks 

We have presented a detailed analysis of the solution of Einstein's equations for 
a static and hyperplane symmetric perfect fluid satisfying the equation of state 
p = r/p, for arbitrary N > 2. 

For r] > — 1, this exact solution clearly show how the attraction of distant 
matter can shrink the spacetime in such a way that it finishes at an empty singular 
boundary. Therefore, in spite of the weakening of gravity with the number of 
dimensions of the space-time [6l[7] , our solution clearly shows that it is still strong 
enough to generate empty singular boundaries for arbitrary N > 2. 

For >)>lwe studied the properties of the geodesies, we found that only vertical 
null geodesies touch the boundary and bounce, and all of them start and finish at 
z = oo; whereas non- vertical null as well as all time-like ones are bounded between 
two planes determined by initial conditions. 

We also showed that the Cauchy problem for the propagation of a massless 
scalar field is well-posed, if we only demand the waves to have finite energy, al- 
though no boundary condition is required. And waves are completely reflected at 
the singularity. 



A Killing vectors and adapted coordinates 

We want to find coordinates adapted to a hyperplane symmetric distribution of matter. That 
is, spacetime must be invariant under n translations and under rotations in n(N — l)/2 hy- 
perplanes. 

More precisely, a n + 2 dimensional spacetime will be said to be n-dimensional Euclidean 
homogenous if it admits the r = n(N + l)/2 parameter group of isometrics n-dimensional 
Euclidean space ISO(n). 

Since the spacetime admits n mutually conmuting independent motions, we can choose 
coordinates x 1 , z, t so that the corresponding Killing vectors are = di (i = 1, . . . , n), and 
so 

( V'i;:::;r) • ^ 

The equations of Killing, 

£, k d k gij + g kj d l £ k + g lk d^ k = , (57) 
corresponding to these vectors reduce to 

fc W =o ( i ' j r=i';;''n z ) ■ («) 
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Hence all the components of the metric tensor depend only on the coordinates t and z and the 
metric is unaltered by the finite transformation 

x i -> x i + a i (for i = 1, . . . , n) . (59) 

We can take for the remaining n(N — l)/2 motions the generators = — x 3 8i (i < 

j and i,j = 1, . . . , n), so 

Taking into account 1581 and (I60D from the equations of Killing I I57II we get 

9jm$i - ffira^ + 9ij<5m - 9iA = 
From the last equation wc readily find 

gu = 9jj and g t j = g it = g iz = 

Furthermore, we can take the bidimensional metric of the Y% spaces x l = constant [i = 1, . . . , ri) 
in the conformal flat form. Hence, the most general metric admitting this group of isometries 
may be written as 

ds 2 = -e™(*& (dt 2 - dz 2 ) + e av 0.*) ((dx 1 ) 2 + ■■■ + (dx n ) 2 ) . (63) 

If, in addition we impose staticity, U and V must be time independent, and the change of 
variable f eF^dz — ► z brings the line element to the form l T3t . 

B Limit cases 

In this appendix, we discuss the limit cases which were not considered section [2] 



(,<7^Vi=C*»)- (61) 

( i j and i,j = 1, , n.) . (62) 



B.l r\ = 1 

By setting e = r? — 1, a = ea and /3 = eb in d 2 5 1) . in the limit e — ► 0, we readily get 

ds 2 = -e 2bu dt 2 + uTT ((dx 1 ) 2 + h (dx N ) 2 ) + — e 2ira u'T d« 2 

4a 

— oo < t < oo, — oo < x 1 < co, . . . , -co < x N < oo, < u < oo . (64) 
Similarly, we get from 1261 1 

p( V ) = p( V ) = ( - N -^ a e -2bu u 4^i (65) 

and from 1281 

„2 e -ibu 

r (N 2 - 1 + 8N(N - l)bu + 16N(N + 2) b 2 u 2 ) . (66) 

6 2 AT2 « 1+ lv v ' 



K : 



Then, since a, b > in this case, density and pressure vanish at u = and u = oo, and the solu- 
tion has an empty singularity at u = 0. This solution is the higher-dimensional generalization 
of that of [16]. 
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B.2 77 = 

It directly follows from i ]2-">[) that, for r) = 0, depending on the sign of a and /3, we get two very 
different solutions 

ds 2 = -n~ 2 3JV+i (1 - /3u) 2 dt 2 +mW+T ((dx 1 ) 2 H h (dx N ) 2 ) + — u"35TI du 2 , 

4a 

i x , < 1//3 if < a, /3 < co 

-co < t < co, — oo < a; < oo, . . . , -co < x < co, < u < lb") 

1 < oo if - oo < a, /3 < . 

In the first case, the solution has all the same properties as solutions in region (III) listed in 
Tabled In the second case, it has all the same ones as those in (II). 

The matter, in this case, is a weird fluid with null density but non-vanishing pressure, and 
H26I) becomes 

N — 1 

p(«) = 0, p(u) = ~ \ ' - — . (68) 

(3N + 1)ttG (1 — pu) 



B.3 77 = -1 

In this case, {25) becomes 



ds 2 = -«-N+t (i _ ^ 2 +7(«tT ((ds 1 ) 2 H h (dx N ) 2 ) + — du 2 , 

V ' U ' y 4a u(l- /3u) 

-oo < t < co, -co < a: 1 < co, . . . , -oo < i™ < oo, < u < 1//3 . (69) 
From (126 ft it follows 

p{u) = -p{u) = -4 - , 70 

which shows that d 69 II turns out to be a vacuum solution with cosmological constant A = 
(jv+i) a < * ms * s * ne higher-dimensional generalization of that found in 1121 . In this case the 
integral in H27H can be computed in closed form, and we get that the depth of this spacetime, 
independently of /3 and N, is 2 ^- . 

On the other hand, we get from (1281) 



16 a 2 (N 2 {N - 1) + 2(N + 2) u 2 /3 2 ) 



K {N + lfu 2 l3 2 ' (71) 

16 a 2 ( 4+JV(JV-2) J 

then, there is a singularity at u = 0, and K = (jv+i) 2 at = 1//3. 



B 4 n - 3N+1 

We now set e = rj(N — 1) + 37V + 1, a = ea and = efe, and by making the coordinate change 
V = 1 — /3it and rescaling t, x' and a, we get the limit e — > of 125t 

ds 2 = _ e 2 (iV-D bv y ^ M 2 + £ -4bv ((^1)2 + . . . + ((faW)2 N + A e - 2 (iV+D 6, U -S^l ^2 

4a 

—co < t < oo, — oo < x 1 < oo, . . . , — oo < x N < oo, < v < oo . (72) 

Similarly, (1261) becomes 

, , 3jV + l , , a(7V-l) 2 2 (N+i)bv ^±1 , m , 
p(v) = p(f ) = — e 2( " +1 > bv v 2N , (73) 
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and (1286 gives 




) ■ (74) 



Therefore, since a, b > 0, density and pressure vanish at v = 0, but diverge at v = oo; and the 
solution has a singularity at v = oo. 
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